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In recent years, there are attempts to describe quantization of planetary distance 
based on time-independent gravitational Schrodinger equation, including Rubcic & 
Rubcic’s method and also Nottale’s Scale Relativity method. Nonetheless, there is 
no solution yet for time-dependent gravitational Schrodinger equation (TDGSE). In 
the present paper, a numerical solution of time-dependent gravitational Schrodinger 
equation is presented, apparently for the first time. This numerical solution leads 
to gravitational Bohr-radius, as expected. In the subsequent section, we also discuss 
plausible extension of this gravitational Schrodinger equation to include the effect 
of phion condensate via Gross-Pitaevskii equation, as described recently by Moffat. 
Alternatively one can consider this condensate from the viewpoint of Bogoliubov- 
deGennes theory, which can be approximated with coupled time-independent 
gravitational Schrodinger equation. Further observation is of course recommended 
in order to refute or verify this proposition. 



1 Introduction 

In the past few years, there have been some hypotheses sug- 
gesting that quantization of planetary distance can be derived 
from a gravitational Schrodinger equation, such as Rubcic 
& Rubcic and also Nottale’s scale relativity method [1, 3], 
Interestingly, the gravitational Bohr radius derived from this 
gravitational Schrodinger equation yields prediction of new 
type of astronomical observation in recent years, i.e. extra- 
solar planets, with unprecedented precision [2]. 

Furthermore, as we discuss in preceding paper [4], using 
similar assumption based on gravitational Bohr radius, one 
could predict new planetoids in the outer orbits of Pluto 
which are apparently in good agreement with recent observa- 
tional finding.. Therefore one could induce from this observ- 
ation that the gravitational Schrodinger equation (and gravi- 
tational Bohr radius) deserves further consideration. 

In the meantime, it is known that all present theories 
discussing gravitational Schrodinger equation only take its 
time-independent limit. Therefore it seems worth to find out 
the solution and implication of time-dependent gravitational 
Schrodinger equation (TDGSE). This is what we will discuss 
in the present paper. 

First we will find out numerical solution of time-inde- 
pendent gravitational Schrodinger equation which shall yield 
gravitational Bohr radius as expected [1, 2, 3]. Then we ex- 
tend our discussion to the problem of time-dependent grav- 
itational Schrodinger equation. 

In the subsequent section, we also discuss plausible ex- 
tension of this gravitational Schrodinger equation to include the 



effect of phion condensate via Gross-Pitaevskii equation, 
as described recently by Moffat [5]. Alternatively one can 
consider this phion condensate model from the viewpoint of 
Bogoliubov-deGennes theory, which can be approximated 
with coupled time-independent gravitational Schrodinger 
equation. To our knowledge this proposition of coupled time- 
independent gravitational Schrodinger equation has never 
been considered before elsewhere. 

Further observation is of course recommended in order 
to verify or refute the propositions outlined herein. 

All numerical computation was performed using Maple. 
Please note that in all conditions considered here, we use 
only gravitational Schrodinger equation as described in Rub- 
cic & Rubcic [3], therefore we neglect the scale relativistic 
effect for clarity. 

2 Numerical solution of time-independent gravitational 
Schrodinger equation and time-dependent gravita- 
tional Schrodinger equation 

First we write down the time-independent gravitational 
Schrodinger radial wave equation in accordance with Rubcic 
& Rubcic [3]: 

d 2 R 2 dR 8Ttm 2 E l 

dr 2 r dr + H 2 , I , 

2 4ir 2 GMm 2 1(1 + 1) 

H R ^ R — 0 . 

r H 2 r 2 

When H, V, E' represents gravitational Planck constant, 
Newtonian potential, and the energy per unit mass of the 
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orbiting body, respectively, and [3]: 




*=* H^r)' 


(2) 


V(r) = - GMm , 


(3) 


E<=° r . 

m 


(4) 


By assuming that R takes the form: 




R = e~ ar 


(5) 



and substituting it into equation (1), and using simplified 
terms only of equation (1), one gets: 

_ ar 2a e~ ar 8irGMm 2 e~ ar 

* = a e e ar + — . (6) 

r r H 2 

After factoring this equation (6) and solving it by equat- 
ing the factor with zero, yields: 



2 UirGMm 2 - H 2 a ) 

RR = 2772 = 0 ■ 

a 2 H 2 



or 



(7) 

( 8 ) 



RR = 4tt GMm 2 - H 2 a = 0 
and solving for a, one gets: 

4ir 2 GMm 2 
a ~ H 2 ' 

Gravitational Bohr radius is defined as inverse of this 
solution of a, then one finds (in accordance with Rubcic & 
Rubcic [3]): 

H 2 



(9) 



ri 4 ir 2 GMm 2 ’ (10) 

and by substituting back equation (2) into (11), one gets [3]: 

27t/\ 2 



r i = 



GM, 



ac 



(11) 



which is equivalent with Nottale’s result [1, 2], especially 
when we introduce the quantization number: r n = r\ n 2 [3], 
For complete Maple session of these all steps, see Appen- 
dix 1. 

Solution of time-dependent gravitational Schrodinger 
equation is more or less similar with the above steps, except 
that we shall take into consideration the right hand side 
of Schrodinger equation and also assuming time dependent 
form of r: 

fl = e-"( f ). (12) 

Therefore the gravitational Schrodinger equation now 
reads: 

d 2 R 2 dR 8irm 2 E ' 

dr 2 + r dr H 2 

2 Aw 2 GMm 2 £(£+1) dR (13) 

ZJ2 ^ 2 ^ ^ j # ’ 

r H z r z at 



or by using Leibniz chain rule, we can rewrite equation 
(15) as: 

dR dr(t ) d 2 R 2 dR 8irm?E' 

dr(t ) dt + dr 2 r dr + H 2 , . 

(14) 

2 47 r 2 GMm 2 _ 1(1+ 1) „ 

H 775 R 2 R = 0 • 

r H z r 2 

The remaining steps are similar with the aforementioned 
procedures for time-independent case, except that now one 
gets an additional term for RR: 

RR 1 = H 3 a^j t r(t)Jr(t) - a 2 r(t)H 2 + 

+ 8-KGMm 2 — 2 H 2 a = 0 . 



At this point one shall assign a value for ^ r ( t ) term, 
because otherwise the equation cannot be solved. We choose 
|r(f) = l for simplicity, then one gets solution for (17): 



a2:= 



a — a, 7T = 7r, m = m, H = H, G = G, M = M , 
t = RootOf(r(_Z)aH 3 -r(_Z)a 2 H 2 + 8Tr 2 GMm 2 -2aH 2 ) 



{a = 0, t — t, m — m, H = H, G = G, M = M, tt = 0} , 
{a = 0, tt — 7T, t — t, m — m, H — H, M = M, G = 0} , 
I tt = tt t — t, m — m, H — H, M — M, a. — H, G — — 

l ’ ’ 4 tt 

{a = a, H — 0, 7T = 7r, t — t, m = m, M — M, G = G} , 
{a = 0, 7r = 7r, t — t, m — m, H — H, G — G, M = 0} , 
{a = 0, 7r = 7r, t = t, H — H, G = G, M = M, m = 0} , 
{a = a, H — 0, 7r = 7r, t — t, m — m, G = G, M = 0} , 
{a = a, H — 0, 7T = 7T, t — t, G = G , M — M, m = 0} 



Therefore one can conclude that there is time-dependent 
modification factor to conventional gravitational Bohr radius 
solution. For complete Maple session of these steps, see 
Appendix 2. 



3 Gross-Pitaevskii effect. Bogoliubov-deGennes appro- 
ximation and coupled time-independent gravitational 
Schrodinger equation 



At this point it seems worthwhile to take into consideration a 
proposition by Moffat, regarding modification of Newtonian 
acceleration law due to phion condensate medium, to include 
Yukawa type potential [5, 6]: 



a(r) = - 



GooM exp (-/i 0 r) 

— 7 — + ^ 7-^(1 +H T )- ( 16 ) 



Therefore equation (1) can be rewritten to become: 
8irm 2 E l 



d?R 2 dR 
dr 2 r dr 



H 2 



-R- 



2 47 t 2 (GM - K exp(— jU0r)(l + fj,^,r))m 2 
+ r If 2 



R - (17) 



1(1+ 1) 



R = 0, 



or by assuming fi = 2 hq = /io r for the exponential term, eq- 
uation (17) can be rewritten as: 
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d?R 

dr 2 



2 dR 
r dr 



8-Km 2 E' 
H 2 



-R- 



2 4t r 2 (GM — ife~ 2/i °(l 
r H 2 



Mo r)) 



R- 



(18) 



Then instead of equation (8), one gets: 

RR" = 87rGMm 2 -2H 2 a-8ir 2 m 2 Ke~^°(l+fj,)=0. (19) 



Solving this equation will yield a modified gravitational 
Bohr radius which includes Yukawa effect: 



H 2 

4ir 2 {GM - Ke- 2 ^)m 2 



( 20 ) 



and the modification factor can be expressed as ratio be- 
tween equation (20) and (11): 



GM 

{GM - Ke- 2 ^°) ’ 



( 21 ) 



for complete Maple session of these steps, see Appendix 3. 

A careful reader may note that this “Yukawa potential 
effect” as shown in equation (21) could be used to explain 
the small discrepancy (around ±8%) between the “observed 
distance” and the computed distance based on gravitational 
Bohr radius [4, 6a]. Nonetheless, in our opinion such an 
interpretation remains an open question, therefore it may be 
worth to explore further. 

There is, however, an alternative way to consider phion 
condensate medium, i.e. by introducing coupled Schrodinger 
equation, which is known as Bogoliubov-deGennes theory 
[7]. This method can be interpreted also as generalisation of 
assumption by Rubcic-Rubcic [3] of subquantum structure 
composed of positive-negative Planck mass. Therefore, 
taking this proposition seriously, then one comes to hypo- 
thesis that there shall be coupled Newtonian potential, in- 
stead of only equation (3). 

To simplify Bogoliubov-deGennes equation, we neglect 
the time-dependent case, therefore the wave equation can be 
written in matrix form [7, p. 4]: 



M [*] =o, 



( 22 ) 



where A] is 2x2 matrix and [SI/] is 2x1 matrix, respecti- 
vely, which can be represented as follows: 

( 8tt GMm 2 e- ar _-ar 2a e - Qr \ 
rTi 2 

2 ae- ar 



M = 



V 



87T GMm e 
rTi 2 



(23) 



and 



M = 



fi r ) 

9{r) 



(24) 



Numerical solution of this matrix differential equation 
can be found in the same way with the previous methods, 
however we leave this problem as an exercise for the readers. 

It is clear here, however, that Bogoliubov-deGennes ap- 
proximation of gravitational Schrodinger equation, taking 
into consideration phion condensate medium will yield non- 
linear effect, because it requires solution of matrix differen- 
tial equation* (22) rather than standard ODE in conventional 
Schrodinger equation. This perhaps may explain complicated 
structure beyond Jovian Planets, such as Kuiper Belt, inner 
and outer Oort Cloud etc. which of course these structure 
cannot be predicted by simple gravitational Schrodinger eq- 
uation [1, 2, 3]. In turn, from the solution of (22) one could 
expect that there are multitude of celestial objects not found 
yet in the Oort Cloud. 

Further observation is also recommended in order to 
verify and explore further this proposition. 

4 Concluding remarks 

In the present paper, a numerical solution of time-dependent 
gravitational Schrodinger equation is presented, apparently 
for the first time. This numerical solution leads to gravita- 
tional Bohr-radius, as expected. 

In the subsequent section, we also discuss plausible ex- 
tension of this gravitational Schrodinger equation to include 
the effect of phion condensate via Gross-Pitaevskii equation, 
as described recently by Moffat. Alternatively one can con- 
sider this condensate from the viewpoint of Bogoliubov- 
deGennes theory, which can be approximated with coupled 
time-independent gravitational Schrodinger equation. 

It is recommended to conduct further observation in order 
to verify and also to explore various implications of our pro- 
positions as described herein. 
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Appendix 2 Time-dependent gravitational Schrodinger equation 

> #Solution of gravitational Schrodinger equation (Rubcic, Fizika 1998); 

> restart; 

> #with time evolution (Hagendom’s paper); 

> S:=r(t); R:=exp(— (alpha*S)); Rl:=exp(— (alpha*r)); 

5 := r(t) 

R := e~ ar 



> D4R:=diff(S,t); D1R:=— alpha*exp(— (alpha*S)); D2R:=— alpha"2* 
exp (— (alpha*S)); D5R:=D1R*D4R; 

°4 R :=|r (t) 

DlR:=-oce~ arW 



DIR : = 



2 



—a 



e ~a :r(t) 



DIR := -ae “ r(t) — r (t) 
dt 



> #Using simplified terms only from equation (A*8) 

> SCHEQ3:= -h*D5R+D2R+DlR*2/S+8V~2*G*M*nT2*R/(S*lT2); 

> XX2:=factor(SCHEQ3); 



XX2 



e -ar(t) ^H 3 a^AIlr(t)-a 2 r(t)H 2 -2H 2 a + 8iT 2 GMm 2 ^ 
r(t)H 2 



> #From standard solution of gravitational Schrodinger equation, we know (Rubcic, 
Fizika 1998): 

> SCHEQ4:=4*pr2*G*M*nT2— !T2*alpha; 



SCHEQ4 := 4ir 2 GMm 2 - H 2 a 



Appendix 1 Time-independent gravitational Schrodinger equation 

> restart; 

> with (linalg); 

> R: = exp(— (alpha*r)); 

R:=e~ ar 

> DlR:=diff(R,r); D2R:=diff(DlR,r); 

DIR := —a e~ a r 
D2R := -a 2 e _Q:r 



> #Therefore time-dependent solution of Schrodinger equation may introduce new 
term to this gravitational Bohr radius. 

> SCHEQ5 : =(XX2*(S *!T2)/(exp ( — (alpha* S)))) — 2* SCHEQ4; 

ODESCHEQ5 := H 3 a dl ^-r(t) - a 2 r(t)H 2 
dt 

> #Then we shall assume for simplicity by assigning value to d[r(t)]/dt: 

> D4R:=1; 

> Then we can solve again SCHEQ5 similar to solution of SCHEQ4: 

> a2 : = solve((lT3 * alpha * (D4R) * S — alpha"2 * S *1T2) +2 * SCHEQ4) ; 



> SCFIEQl:=D2R+DlR*2/r+8V'2*m*E*R/h'2+8*pr2*G*M*nr2*R/(r*iT2)— 
l*(l+l)*R/r~2=0; 

> XX 1 : = factor (SCHEQ 1 ); 

> #Using simplified terms only from equation (A*8, of Rubcic & Rubcic, 1998) 

> ODESCHEQ:=D2R+DlR*2/r+8*pr2*G*M*m''2*R/(r*h"2)=0; 



ODESCHEQ := a 2 e~ ar 



2a e~ a r 8ir 2 GMm 2 e- a r 
r + rH 2 



= 0 



> XX2: =factor (SCHEQ2); 



XX2 



e~ ar ( a 2 rH 2 - 2 H 2 a + 8tt 2 GMm 2 ) 
rH 2 



= 0 



a — a, 7r = 7 r, m = m, H — H, G — G, M = M, 
t = RootO f(r(_Z)aH 3 -r(_Z)a 2 H 2 +8Tr 2 GMm 2 -2aH 2 

{a = 0, t — t, m — m, H = H , G — G, M = M, 7r = 0} , 

{a = 0, 7r = 7T, t — t, m — m , H = H, M = M, G = 0} , 

S tt = tt, t = t, m — m, H = H, M = M, a = H, G = 3 9 \ , 

\ 4tt 2 Mm 2 J ’ 

{a = a, H = 0, 7r = 7r, t — t, m — m, M — M, G = G}, 

{a = 0, 7r = 7r, t — t, m — m, H = H , G — G, M = 0} , 

{a = 0, 7r = 7T, t — t, H — H, G = G, M — M, m = 0} , 

{a = a, H — 0, 7r = 7r, t — t, m = m, G — G, M = 0} , 

{a = a, H — 0, 7T = 7r, t — t, G — G, M — M, m = 0} 




> RR:= solve (XX2,r); 



RR := - 



2(4t r 2 GMm 2 - H 2 a) 
a 2 H 2 



> #Then solving for RR=0, yields: 

> SCHEQ3: = 4*pr2*G*M*m A 2-r2*alpha=0; 

SCHEQ3 := 4ir 2 GMm 2 - H 2 a = 0 

> a:= solve (SCHEQ3, alpha); 

4ir 2 GMm 2 

a := 

H 2 

> #Gravitational Bohr radius is defined as inverse of alpha: 

> gravBohrradius: = l/a; 

_ H 2 

rgravBohr 



> #Therefore one could expect that there is time-dependent change of gravitational 
Bohr radius. 

Appendix 3 Time-independent gravitational Schrodinger equation 
with Yukawa potential [5] 

> #Extension of gravitational Schrodinger equation (Rubcic, Fizika 1998); 

> restart; 

> #departure from Newton potential; 

> R : = exp ( — (alpha*r)) ; 

R-e~ ar 

> DlR:=diff(R,r); D2R:=diff(DlR,r); 



DIR := - ae ar 
D2R := -a 2 e~ ar 
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> SCHEQ2:=D2R+DlR*2/r+8V"2*(G*M— K*exp(— 2*mu)*(l+mu*r))*nf2*R/ 
(r*h"2)=0; 

2 2a e~ a r 

ODESCHEQ := a 2 e ar + 

r 

8-k 2 (GM - Ke~ 2lJ -(l + ^r))m 2 e~ ar _ n 

+ Vh 2 “ ° 

> XX2:=factor(SCHEQ2); 

> RRl:=solve(XX2,r); 

2 f-H 2 a + 4n 2 GMm 2 - 4n 2 m 2 ife -2 ") 

—a 2 H 2 + 8ir 2 m 2 Ke~ 2 ^ 



> #from standard gravitational Schrodinger equation we know: 

> SCHEQ3:=4*pr2*G*M*m~2— h A 2*alpha=0; 

> a:=solve(SCHEQ3,alpha); 

> #Gravitational Bohr radius is defined as inverse of alpha: 

> gravBohrradius: = l/a; 

H 2 

r gray Bohr ■- ^2 GMm 2 

> #Therefore we conclude that the new terms of RR shall yield new terms (YY) into 
this gravitational Bohr radius: 

> PI:= (RR*(alpha~2*h"2)-(-8 V~2*G*M*m A 2+2*lY2*alpha)); 

> #This new term induced by pion condensation via Gross-Pitaevskii equation may 
be observed in the form of long-range potential effect, (see Moffat J., arXiv: astro- 
ph/0602607, 2006; also Smarandache F. and Christianto V. Progress in Physics, v. 2, 
2006, & v. 1, 2007, www.ptep-online.com) 

> #We can also solve directly: 

> SCHEQ5 : = RR* (alpha"2 *!T2)/2 ; 



SCHEQ5 : = 



a 2 H 2 (-H 2 a + 4-n 2 GMm 2 - 4ir 2 m 2 Ke~ 2lJ -) 
-a 2 H 2 + 8 ? r 2 m 2 Ke- 2 ^ 



> al:=solve(SCHEQ5, alpha); 



a 1 := 0,0, 



4ir 2 m 2 (GM - 
H 2 



> #Then one finds modified gravitational Bohr radius in the form: 

> modifgravBohrradius:= l/(4*pf2*(G*M— K*exp (-2*mu))*nY2/h"2); 

H 2 

r modified, gray Bohr ■- 4 ^2 m 2 (GM _ Ke -2» j 



> #This modification can be expressed in chi-factor: 

> chi: =modifgravBohrradius/gravBohrradius; 

GM 

X GM - Ke- 2 * 
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